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Abstract 

We give a brief account of recent developments in the theoretical understanding of the interac- 
tion between electric currents and inhomogeneous ferromagnetic order parameters. We start by 
discussing the physical origin of the spin torques responsible for this interaction and construct a phe- 
nomenological description. We then consider the electric current-induced ferromagnetic instability 
and domain-wall motion. Finally, we present a microscopic justification of the phenomenological 
description of current-driven magnetization dynamics, with particular emphasis on the dissipative 
terms, the so-called Gilbert damping a and the (3 component of the adiabatic current-driven torque. 
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I. INTRODUCTION 



Ferromagnetism is a correlated state in which, at sufficiently low temperatures, the elec- 
trons align their spins in order to reduce the exchange energy. Below the Curie temperature 
Tc, the free energy F[M] then attains its minimum at a finite magnetization M 7^ with 
an arbitrary direction, thus spontaneously breaking the spin-rotational symmetry. Crystal 
anisotropics that are caused by the spin-orbit interaction and shape anisotropics governed 
by the magnetostatic dipolar interaction pin the equilibrium magnetization direction to a 
certain plane or axis. At low temperatures, T <C T^, fiuctuations of the magnitude of the 
magnetization around the saturation value (the so-called Stoner excitations) become 
energetically unfavorable. The remaining low-energy long-wavelength excitations are spin 
waves (or magnons, which, in technical terms, are viewed as Goldstone modes that restore 
the broken symmetry). These are slowly varying modulations of the magnetization direction 
in space and time. 

A phenomenological description of the slow collective magnetization dynamics without 
dissipation proceeds from the free energy F[M(r)] as a functional of the inhomogeneous 
(and instantaneous) magnetic configuration M(r) pL]- The equation of motion 

(9M(r,t) , 5F\M] 

^^ = 7M(r,t)x — ^ — ^ (1) 

preserves the total free energy of the system, since the rate of change of the magneti- 
zation is perpendicular to the "gradient" of the free energy. The functional derivative 
of the free energy with respect to the local magnetization is called the effective field: 
Heff(r, t) = — (5F[M]/5M. [In the following, we use the abbreviations dt = d/dt for the 
partial derivative in time and (9m = for the functional derivative with respect to M.] 

In the presence of only an externally applied magnetic field H, F[M] = — J (i^rM(r) ■ H(r), 
so 7 is identified as the effective gyromagnetic ratio. In general, Hefj includes the crys- 
tal anisotropy due to spin-orbit interactions, modulation of the exchange energy due to 
magnetization gradients, and demagnetization fields due to dipole-dipole interactions. The 
^andau-Lifshitz equation ([1]) qualitatively describes many ferromagnetic resonance (FMR) 
2I, [3] and Brillouin light scattering j4| experiments. With M = M^m, where m is the 
magnetic direction unit vector, we may rewrite Eq. ([1]) as 

dMr, t) = -7m(r, t) x Heff(r, t) . (2) 



By definition, the effective field on the right-hand side of Eq. ([2]) is determined by the 
instantaneous magnetic configuration. This can be true only if the motion is so slow that 
all relevant microscopic degrees of freedom manage to immediately readjust themselves to 
the varying magnetization. If this is not the case, the effective field acquires a finite time 
lag that to lowest order in frequency can be schematically expanded as 



Heff ^ -dMF[M{T, t-r)]^ Heff - r{dtM ■ 5M)Heff . 



(3) 



where r is a characteristic delay time. This dynamic correction to the instantaneous effective 
field, 5Hefj, leads to a new term in the equation of motion oc m x 5Hcfj. Although in general 
nonlocal and anisotropic j^, it makes sense to identify first the simplest, i.e., local and 
isotropic contribution. We can then construct two new terms out of the vectors m and 
dtm. The first one, oc m x dtui, is dissipative, meaning that it is odd under time reversal 
(i.e., under the transformation t — > — t, H — H, and m — m), and thus violates the 
time-reversal symmetry of the Landau-Lifshitz equation ([2]). This argument leads to the 
Landau-Lifshitz-Gilbert (LLG) equation [g, l7|: 



9tm(r, t) = — 7m(r, t) x Hefr(r,t) + am(r, t) x dtm{r,t) 



(4) 



introducing the phenomenological Gilbert damping constant a. The second local and 
isotropic term linear in dtin and perpendicular to m, that can be composed out of m 
and dtm is proportional to dtm and can be combined with the left-hand side of the LLG 
equation. In principle, any physical process that contributes to the Gilbert damping can 
thus also renormalize the gyromagnetic ratio 7. The latter should therefore be interpreted 
as an effective parameter in the equation of motion (jl]). The second law of thermodynamics 
requires that 07 > 0, which guarantees that the dissipation of energy P oc Hcs • dtm > 
(assuming the magnetization dynamics are slow and isolated from any external sinks of 
entropy). Since the implications of small modifications of the gyromagnetic ratio are mi- 
nor, we will be mainly concerned with the Gilbert damping constant a. Furthermore, we 
will focus on dissipative effects due to spin dephasing by magnetic or spin-orbit impurities 

notiiig that many other Gilbert damping mechanisms have been pro- 




posed in the past 



20 



21| . According to the fluctuation-dissipation 



theorem, the dissipation, whatever its microscopic origin is, must be accompanied by a 
stochastic contribution h(r, t) to the effective field. Assuming Gaussian statistics with a 



white noise correlator [22|, which is vahd in the classical limit with characteristic frequen- 
cies that are sufficiently small compared to thermal energies: 



a 



{h,{r,t)h,{T',t')) = 2kBT—-6,,6{r - r')6{t - t') . 



(5) 



Eqs. ([I])-© form the standard phenomenological basis for understanding dynamics of fer- 
romagnets , in the absence of an applied current or voltage bias. 



II. CURRENT-DRIVEN MAGNETIZATION DYNAMICS 
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n order to understand recent experiments on current-biased magnet ic multilayers 
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modified 



27 



28 



29 



30 



311] and nanowires 32, 



33 
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39|, Eq. 



has to be 



40, |4l|, |42] . The leading correction has to take into account the finite divergence 



of the spin-current density in conducting ferromagnets, with magnetization texture that has 
to be brought into compliance with the conservation of angular momentum. We have to 
introduce a new term 

SO^i^^iltorque = V • ji (6) 

in the presence of a current density jj for spin-z component, where Sq is the total equilibrium 
spin density along — m (the minus sign takes into account that electron spin and magnetic 
moment point in opposite directions) . By adding this term to the right-hand side of Eq. (jlj) 
as a contribution to dtrrii, we assume that the angular momentum lost in the spin current 
is fully added to the magnetization. This is called the spin-transfer torque 42|, |43|, 14^, |45 |. 
The simplest approximation for the spin-current density in the bulk of isotropic ferromagnets 
VjiTLi, where P is a material-dependent constant that converts charge- 



43 



45 



m is j^ 



current density j into spin-current density. The underlying assumption here is that spins are 
carried by the electric current such that the spin-polarization axis adiabatically follows the 
local magnetization direction. This is the case for a large exchange field that varies slowly 
in space. This condition is fulfilled very well in transition-metal ferromagnets. The spin 
conversion factor 

V = 1^1^ (7) 

characterizes the polarization of the spin-dependent conductivity (s =t or s =|) with | 
chosen along — m. Hence 



V 

dtjn = —7m X Hefj H ( j ■ V)m . 

So 



(8) 



where we took into account local charge neutrality by V ■ j = 0. 

The phenomenological Eq. ([8]) is "derived" without taking into account spin relaxation 
processes. Its inclusion requires some care since both Gilbert damping and spin-transfer 



torque are nontrivially affected 



47 



48 



49 



50 



51 



52| . Spin relaxation is generated by 



impurities with potentials that do not commute with the spin density operator, such as 
a quenched random magnetic field or spin-orbit interaction associated with randomly dis- 



tributed non-magnetic impurities 



13 



47|, 



49, 



50 



511 ] ■ In the absence of an applied current 



ections with potentials that mix the spin channels contribute to the Gilbert con- 



si, y, 0, Q: 



50l ]. It is instructive to interpret the right-hand side of the equation 



J, imperi 
stant a 

of motion ([H]) as an analytic expansion of driving and damping torques in V and dt. The 
LLG equation (jl]) corresponds to the most general (local and isotropic) expression for the 



damping to the zeroth order in V and first order in dt- We will not be concerned with 
higher order terms in 9t, since the characteristic frequencies of magnetization dynamics are 
typically small on the scale of the relevant microscopic energies, at least in metallic systems. 
The contribution to the effective field due to a finite magnetic stiffness jl| is proportional 
to V^. In the presence of inversion symmetry, the terms proportional to V cannot appear 
without applied electric currents. In the following, we focus on the current- driven terms lin- 
ear in V, assuming that the spatial variations in the magnetization direction are sufficiently 
smooth to rule out higher-order contributions. The dynamics of isotropic spin-rotationally 
invariant ferroma gne ts can then in general be described by the phenomenological equation 



of motion [471, |49|, [53 1 



V 

-7m X Hcff + am X dtva. H (1 

So 



/5mx) (j • V)m , 



(9) 



in which a and /3 characterize those terms that break time-reversal symmetry. Both arise 



47 



49 



50|. Even though in 



47 



practice 



49 



53|, 



naturally in the presence of spin-dependent impurities 

/? <^ 1, it gives an important correction to the current-driven spin-transfer torque 
as discussed in more detail below. For the special case of a = /3: Eq. ([9]) can then be rewritten 
(after multiplying it by 1 + amx on the left) as 



V 

dtui = — 7*m X Heff + a;7*m x Heg x m H (j ■ V)m . 

So 



(10) 



where 7* = 7/(1 + a^). The dissipative term proportional to m x Hefj x m is called the 
Landau-Lifshitz damping. Eq. ([9]) cannot be transformed into equation (|T0|) if a 7^ /3 [in 



which case Eq. (fTOj) necessarily retains a j3 term]. A special feature of Eq. (fTO!) appears when 
Hefr(m) is time independent and translationally invariant. A general solution m(r,t) in the 
absence of an electric current j = (such as a static domain wall or a spin wave) can be 
used to construct a solution 

m(r,t) = m(r + Pjt/so,t) (11) 

of Eq. ( |TOl) for an arbitrary uniform j. This unique feature of the solutions of Eq. ([9]) only 
arises for a = /3. 



Interestingly, the argument above has been turned around by Barnes and Maekawa 54| . 
who find that Galilean invariance of a system would dictate a = (3. Galilean invariance 
requires the existence of solutions of the form m(r — vt), where m(r) is an arbitrary static 
solution (say, a domain wall) and v is an arbitrary velocity. As explained above, this is 
only possible when a = (3. However, the general validity of the Galilean invariance assump- 
tion for the current- carrying state needs to be discussed in more detail from a microscopic 



point of view. The Galilean invariance argument [5J] implies that the bias-induced electron 
drift exactly corresponds to the domain-wall velocity, since otherwise electron motion would 
persist in the frame that moves with the domain wall. Referring to Eq. ffTTl) . we must, 
therefore, identify v = —V]/sq with the average electron drift velocity in the presence of the 
current j. We argue in the following that this is indeed true in certain special limits, but 
is not generic, however. In the itinerant Stoner model for ferromagnets, the spin-dependent 
Drude conductivity reads cTg oc ^sT^, where and are the densities and scattering times 
of spin s, respectively. When there is no asymmetry between the scattering times, = T|, 
V indeed equals the electron drift velocity and Galilean invariance is effectively fulfilled. In 
general, however, since the spin dependence of wave functions and densities of states for 
the electrons at the Fermi energy lead to different scattering cross sections in conducting 
ferromagnets, the equality between — Pj/so and the average drift velocity disappears. In 
the simplest model of perturbative white-noise impurity potentials, for example, l/r^ oc 
where z/^ is the spin-dependent density of states. Assuming parabolic free-electron bands 
and weak ferromagnets, in which the ferromagnetic exchange splitting is much less than the 
Fermi energy, the domain- wall velocity v = — Pj/so actually becomes 2/3 of the average 
drift velocity v. 

So n^T^ + n^Ti n^-ni {n^ - ni){v^ + u^) /2 3 

6 



1/3 

Here, is the spin-s electron drift velocity, and we used the relation z/^ oc nj , which 
is valid in three dimensions. Clearly, the potential disorder breaks Galilean invariance. 
An identity of a and (3 can therefore not be deduced from general symmetry principles. 
Furthermore, spin-orbit interaction or magnetic disorder that strongly affect the values of a 
and (3 (see below) also break Galilean invariance at the level of the microscopic Hamiltonian. 
Nevertheless, for itinerant ferromagnets we show below that a ~ /3 (where by ~ we mean 
"of the order"), with a j3 in the simplest model of weak and isotropic spin-dephasing 



impurities 



491], which implies that deviations from translational invariance are not very 



important in metallic ferromagnets, such as transition metals and their alloys, in which 
the Stoner model is applicable. Very recently, two independent groups measured a ^ (3 vn 
permalloy nanowires jsa. Is?]. 

Let us also consider the s — d model of ferromagnetism. When, as is usually done, the d- 
orbital lattice is assumed spatially locked, Galilean invariance is broken even in the absence 
of disorder. In this case, the ratio a/ (3 deviates strongly from unity, although it remains to 
be relatively insensitive to the strength of spin-dependent impurities. In other words, a and 
13 scale similarly with the strength of spin-dephasing processes, and their ratio appears to 
be determined mainly by band-structure effects and the nature (rather than the strength) 



of the disorder 



49 



50|]. A predictive material-dependent theory of magnetization damping 



and current-induced domain wall motion that transcends the toy models mentioned above 
is beyond the scope of our paper. 

The form of the equation of motion ( fTOl) for the special case a = (3 has also triggered 
the suggestion 551] that the Landau-Lifshitz form of damping, (X. m X HcfT X m, is more 
natural than the Gilbert form, oc m x 9jm. In our opinion, however, such a distinction 
is purely semantic. Both forms are odd under time reversal, and one can easily imagine 
simple models in which either form arises more naturally than the other: For example, a 
Bloch-like T2 relaxation added to the Stoner model naturally leads to a Landau-Lifshitz 
form of damping j49| , whereas the dynamic interface spin pumping {sgI, [s^] very generally 
obeys the Gilbert damping form. Moreover, mathematically, both equations are identical 
(in the absence of any additional torques), since we have shown above that the Landau- 
Lifshitz form of damping follows from the Gilbert one simply by multiplying both sides of 
the LLG form by 1 + amx from the left (and vice versa by 1 — amx). Only at the special 
point a = /9, the Landau-Lifshitz form ffTOl) does not involve the "/3 term." In that limit, 
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it may be a more transparent expression for the equation of motion. On the other hand, 
we noted above that in general a ^ (3 and the ratio a/ (3 depends on material and sample. 
The current-driven dynamics of domain walls and other spatially nonuniform magnetization 
distributions turn out to be very sensitive to small deviations of a/ (3 from unity, which 
strongly reduces any advantage a Landau-Lifshitz damping formulation might have over the 
Gilbert phenomenology. In general, we therefore prefer to use the Gilbert phenomenology. 

Under time reversal, the electric current as well as the magnetization vector change sign 
and the adiabatic current-induced torque is symmetric, thus nondissipative. The Ohmic 
dissipation generated by this current does not depend on the magnetization texture in this 
limit and is intentionally disregarded. Saslow 58|] prefers to discuss a torque driven by 
voltage rather current, which, after inserting Ohm's law for the current, becomes odd under 
time reversal and thus appears dissipative (see Ref. 59] for another discussion of this point). 
Obviously, the f3 correction torque is odd for the current-biased and even for the voltage- 
biased configurations. The current-bias picture appears to be more natural, since it reflects 
the absence of additional dissipation by the magnetization texture in the adiabatic limit as 
well as the close relation between the f3 correction and the Gilbert dissipation. 



III. CURRENT-DRIVEN INSTABILITY OF FERROMAGNETISM 

Let us now pursue some special aspects of the solutions of the phenomenological Eq. IQ, 
highlighting the role of various parameters, before we discuss the microscopic derivation 
of the magnetization dynamics in Sec. |Vl It is interesting, for example, to investigate 



the 



43 



possibility to destabilize a single-domain ferromagnet by sufficiently large spin torques 



461]. We consider a homogeneous ferromagnet with an easy-axis anisotropy along the x 
axis, characterized by the anisotropy constant K, and an easy-plane anisotropy in the xy 
plane, with the anisotropy constant K±, see Fig. [H Typically, the anisotropies originate 
from the demagnetization fields: For a ferromagnetic wire, for example, the magnetostatic 
energy is lowest when the magnetization is in the wire direction, so that there are no stray 
field lines outside the ferromagnet. The effective field governing magnetization dynamics is 
then given by 

Heff = {H + Km^)x - K±m^z + AV^m , (13) 
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where we also included an applied field H along the x axis and the exchange coupling 
parametrized by the stiffness constant A. A, H, K, K± > 0. We then look for spin-wave 
solutions of the form 

m = X + ue^^'i '-^*) , (14) 

plugging it into Eq. ([9]) with the effective field (fT3i) in the presence of a constant current 
density j, and linearizing it with respect to small deviations u. When a = and j = 0, we 
recover the usual spin- wave dispersion (Kittel formula): 

LOoiq) = lV{H + K + Aq^){H + K + K^ + Aq^) . (15) 

A finite 07 > results in a negative lmuj{q), as required by the stability of the ferromagnetic 
state. A sufficiently large electric current may, however, reverse the sign of Imu;(g) for certain 
wave vectors q, signaling the onset of an instability. The critical value of the current for the 
instability corresponds to the condition lmuj{q) = 0. Straightforward manipulations based 
on Eqs. and (JT3l) show that this condition is satisfied when 

^ ^-^Vq-je) = ±^o(g). (16) 



So V 

which leads to a critical current density 

jcQ is the lowest current satisfying equation ('P/so)(q ■ jco) = ^^o{q) for some q, where the 
left-hand side can be loosely interpreted as the current-induced Doppler shift to the nat- 



ural frequency given by the right-hand side 



461]. According to Eq. f|T7j) . a current-driven 



instability is absent when a = (3. This conclusion is in line with the arguments leading to 
Eq. ffTTl) : For the special case of a = /3, a spin- wave solution in the presence of a finite 
current density j would acquire a frequency boost proportional to q • j, but with a stable 
amplitude. Note that in general the onset of the current-driven ferromagnetic instability 
is significantly modified by the existence of (3 even with /3 <^ 1, provided that the ratio 
P/a is appreciable. In fact, a is typically measured to be ~ 0.001 — 0.01, and the existing 



microscopic theories 



50| predict P to be not too different from a. 
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FIG. 1: Transverse head-to- head (Neel) domain wall parallel to the y axis in the easy xy plane. 
The uniform magnetization has two stable solutions m = ibx along the easy axis x, which is 
characterized by the anisotropy constant K. These are approached far away from the domain wall: 
m ibx at x = Too, respectively. In equilibrium, the magnetization direction m is forced into the 
xy plane by the easy-plane anisotropy parametrized by K^. A weak magnetic field H or electric 
current j applied along the x axis can induce a slow domain- wall drift along the x axis, during 
which the magnetization close to the domain wall is tilted slightly out of the xy plane. At larger H 
or j (above the so-called Walker threshold), the magnetization is significantly pushed out of the xy 
plane and undergoes precessional motion during the drift. In the moving frame, the magnetization 
profile may remain still close to the equilibrium one. 

IV. CURRENT-DRIVEN DOMAIN- WALL MOTION 



Even more interesting phenomena are associated with the effect of the applied electric 
current on a stationary domain-wall. In particular, we wish to discuss how the spin torques 
move and distort a domain wall. These questions date back about three decades [4 
alt houg h orily relatively recently they sparked an intense activity by several groups 



5o[ 

32, 



51 
331 



53 
34 



54 
35 



601. 



36 



491, 



38 



lis is motivated by the growing number of intriguing experiments 



39 1 as well as the promise of practical potential, such as in the so- 



called racetrack memory 6l| or magnetic logics 62]. Current-induced domain- wall motion 
is a central topic of the present review. 

For not too strong driving currents and in the absence of any significant transverse dy- 
namics, one can make progress analytically by using the one- dimensional Walker ansatz. 



which was first employed in studies of magnetic-field driven domain-wall dynamics 



63| 
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This approach has proven useful in the present context as well [60|, ISJ, |65|] . The key idea 
is to approximately capture the potentially complex domain-wall motion by few parameters 
describing the displacement of its center and a net distortion of the domain-wall structure. 
In a quasi-one-dimensional set-up, such as a narrow magnetic wire, the domain wall is con- 
strained to move along a certain axis, whereas the transverse dynamics are suppressed. This 
regime is relevant for a number of existing experiments, although it should be pointed out 
that the common vortex-type domain walls do not necessarily fall into this category. Let us 
consider an idealized situation with an effective field (fT3|l and an equilibrium domain wall 
magnetization in the xy plane. The magnetization prefers to be collinear with the x axis 
due to the easy-axis anisotropy K. A transverse head-to-head domain wall parallel to the y 
axis corresponds to a magnetization direction that smoothly rotates in the xy plane between 
X at a; ^ — oo and — x at s ^ cxd, as sketched in Fig. [H 

The collective domain- wall dynamics can be described by the center position X{t) and 
an out-of-plane tilting angle ^{t). [F or a more technically-interested reader, we note that 



in the effective treatment of Ref. 60|], these variables are canonically conjugate.] There is 
also a width distortion, but that is usually considered less important. When H < K, the 
two uniform stable states are m = ±x. When if = 0, a static transverse head-to- head 
domain-wall solution centered at x = is given by 

^{x)^0, lntan^ = |,, (18) 

where position-dependent angles cp and 6 parametrize the magnetic configuration: 

m = {rrix, my, irtz) = (cos 6, sin 6 cos if, sin 6 sin ip) . (19) 

W = a/ A/K is the wall width, which is governed by the interplay between the stiffness A 
that tends to smooth the wall extent and the easy-axis anisotropy K that tends to sharpen 
the wall. 

The external magnetic field H or the current density j along the x axis disturb the static 
solution ( |T8i) . distorting the domain- wall structure and displacing its position. At weak field 



and current biases, magnetic dynamics can be captured by the Walker ansatz 63|, l6J] : 



. X / X , 0(r,t) X - X(t) , , 

r,t = $ t , Intan^^ = ^ — ^. (20) 

Here, it is assumed that the driving perturbations {H and j) are not too strong, such that 
the wall preserves its shape, except for a small change of its width W{t) and a uniform 

11 



out-of-plane tilt angle X{t) parametrizes the net displacement of the wall along the 

X axis. Note that although ip is assumed to be spatially uniform, it has an effect on the 
magnetization direction only when m 7^ ±x, i.e., only near the wall center. A more detailed 



discussion concerning the range of validity of this approximation can be found in Ref. [63 1. 
Inserting the ansatz (1201) into the equation of motion with j = jx (since the other current 
directions do not couple to the wall), and using Eq. (fT3|) for the effective field, one finds 



53, 



6^ 



• aX pPj 

$ + ^ = 7^ - — , 
W soW 

X ■ 7if_Lsin2<l> Vj 
— a;$ 



W 2 soW 

W = J ^. (21) 

\ K + K^ sin^ $ ^ ' 

It is easy to verify that the static solution (ITS]) is consistent with these equations when H = 
and j = 0. Two different dynamic regimes can be distinguished based on Eqs. (|^ : When 
the driving forces are weak, a slightly distorted wall moves at a constant speed, X = const, 
and constant tilt angle, $ = (assuming constant H and j). The corresponding Walker 
ansatz (!20|) then actually provides the exact solution, which is approached at long times 
after the constant driving field and/or current are switched on 63|, l6J]. Beyond certain 
critical values of H ot j, called Walker thresholds, however, no solution with constant angle 
$ and constant velocity X exist. Both undergo periodic oscillations in time, albeit with a 
finite average drift velocity {X) 7^ 0. In the spacial case of a = /?, Eqs. f l2T]) are exact at 
arbitrary dc currents when H = 0: According to Eq. (ITT]) , the static domain- wall solution 
(of an arbitrary domain-wall shape) then simply moves with velocity —VJ/sq without any 
distortions. When (3 a, the Walker threshold current diverges when (3 approaches a, 
reminiscent of the critical current ( |T71) discussed in the previous section. 

For subthreshold fields and currents with $(t) — » const as t — > 00, the steady state 
terminal velocity is given by 47| 

• ^HW - /3Vj/so 
V = X[t 00) = . (22) 

a 

In particular, when j = 0, the wall depicted in Fig. [1] moves along the direction of the 
applied magnetic field H in order to decrease the free energy 63|]. Let us in the following 



12 



focus on the current- driven dynamics with H = 0. At a finite but small j, the wall is slightly 
compressed according to 

^ W 2^^AK^ a) ' ^^^^ 



where W = A/K is the equilibrium width. When a = (3, the domain-wall velocity 
V —Vj/sQ. In this case, if we consider the electron spins following the magnetization 
direction from ±m to =Fm on traversing the domain wall with current density j, the entire 
angular momentum change is transferred to the domain- wall displacement. In this sense, 
the ratio (3 /a can be loosely interpreted as a spin-transfer efficiency from the current density 
to the domain-wall motion. 

Only when a = (3, the rigidly moving domain-wall solution is exact at arbitrary current 
densities, leading to an infinite Walker threshold current. The latter becomes finite and 
decreases with /3 < a, approaching a finite value jto at /5 = 6(|, see Fig. [2l In the absence 
of a strong disorder pinning centers, as assumed so far, jto oc K± (which is also the case 



with the Walker threshold field in the absence of an applied current 631]), with an average 
velocity that slightly above the threshold reads 



(X) (X - 3l . (24) 

See the /3 = curve in Fig. [21 At finite /3, the depinning current is determined by the pinning 
fields, which should be included into the effective field ( !T3|) . The domain- wall velocity at 



54| to grow linearly with 



currents slightly above the depinning current is predicted in Ref. 
3- 

So far in our discussion, we have completely disregarded the random noise contribution 
to the magnetization dynamics. As noted above, see Eq. ([5]), thermal fiuctuations are 
ubiquitous in dissipative systems. Below the (zero-temperature) depinning currents, apphed 
currents can drive the domain wall with finite average velocity (v) only by thermal activation. 
The question how ln(f ) scales with the current at low temperatures and weak currents is 
of fundamental interest beyond the field of magnetism. Experiments on thermally-activated 



domain- wall motion in magnetic semiconductors 33|, l39|] reveal a "creep" regime 66|], in 
which the effective thermal-activation barrier diverges at low current density j, so that 
\n.{v) scales as const — with an exponent yU ~ 1/3. This is inconsistent with the 
theory based on the Walker ansatz for rigid domain- wall motion [gT], which yields a simple 
linear scaling of the effective activation barrier and ln(w) oc const + j. A refinement of the 

13 



0.5 1 1.5 2 2.5 3 



FIG. 2: Average current-driven domain-wall velocity v numerically calculated using the Walker 
ansatz [Eqs. ()21|) ] in Ref. [53!]- Here, the domain- wall width has been approximated by its equilib- 
rium value, W ~ W , assuming K± <C K. The curves are very similar to the full micromagnetic 
simulations 53|. u = —Vj/sq has the units of velocity (proportional to electron drift velocity) and 



^K^C,/2 is its value for j = jtQ. The length C, ~ W ^ if we assume i^_L ^ K (as was done 



in this calculation), while C, ~ ^J2A/K_\_ in the opposite limit, K_\_ ^ which is relevant for a 
thin-film with large demagnetization anisotropy K±^ = 4ttMs (in which case C is called exchange 



length) 



64( 1 ■ a = 0.02, and we refer to Ref. 



53] for the remaining details. 



Walker-ansatz treatment 



68 



cannot explain the experiments either. A scaling theory of 



creep motion close to the critical temperature 



measurements by Yamanouchi et al. 



39| does offer a qualitative agreement with 



39j. However, the intrinsic spin-orbit coupling in p- 



doped (Ga,Mn)As 



see, e.g., Refs. 



69 



eads to current-driven effects beyond the standard spin-transfer theories, 



TOj, which needs to be understood better in the present context. 



Even at zero temperature, there are stochastic spin-torque sources in the presence of 
an applied current, which stem from the discreteness of the angular momentum carried by 
electron spins, in analogy with the telegraph-like shot noise of electric current carried by 
discrete particles. A theoretical study of the combined thermal and shot-noise contributions 
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to the stochastic torques for inhomogeneous magnetic configurations 7l| did not yet explore 
consequences for the domain-wall dynamics. For example, it is not known whether shot noise 
assists the current-driven domain-wall depinning at low temperatures. Questions along these 
lines pose challenging problems for future research. 

Effects beyond the theory discussed above are generated by nonadiabatic spin torques, 
which lead to higher-order in dt and V terms in the equation of motion ([9]) . It is in p rinciple 



13 



49, 



50| 



possible to extend linear-response diagrammatic Green's function calculation 
by systematically calculating higher-order terms as an expansion in the small parameters, 
i.e., spin- wave frequency and momentum [t^. A dynamic correction to the spin torque in 



Eq. 



has been found in Refs. 



49|,l72|, which comes down to replacing (3 (3 + n{h/ A^c)dt, 



where is the ferromagnetic exchange splitting and n = 1{2) for the Stoner (s — d) 



model 



49|, 



72| . Since this term scales like dt'V, it is symmetric under time reversal and 



therefore nondissipative. Although this dynamic correction is rather small at the typical 
FMR frequencies huj <^ Axe, it can cause significant effects at large currents [72|]. Starting 



from an inhomogeneous equilibrium configuration 



72 



73 



74 



75| , such as a magnetic spiral or 



a domain wall, one can capture nonadiabatic terms in the equation of motion that vanish 



n m 



13 



49 



5G, 



5l|. 



linear response with respect to the uniform magnetization considered in Refs 

For strongly-inhomogeneous magnetic structures, perturbative expansions around a uni 
form magnetic state fail. For example, for sharp domain walls, the effective equations (1211 



describing wall dynamics and displacement acquire a new term, which can be understood 



cLS cL 

wall 



brce transferred by electrons reflected at the potential barrier caused by the domain 



60 



761]. Electron reflection at a domain wall increases the resistance. Adiabaticity 



implies a vanishing intrinsic domain-wall resistance (see, however, Ref. 69|] for a model with 
strong intrinsic spin-orbit coupling). The force term, becomes important only for abrupt 
walls with width W ~ Axe = fivp/A^c- Such nonadiabatic effects are not expected to be 
strong in metallic ferromagnets, where typically W ^ Axe ~ Ai? (the Fermi wavelength). 
Dilute magnetic semiconductors [such as (Ga,Mn)As] are a different class of materials with 



longer Axe and a strong spin-orbit coupling 



70|. 



In metallic systems, effects of the spin-torque in the most relevant regime of slow dy- 
namics, fkv <ti Axe, with smooth walls, W ^ Axe, and at moderate applied currents is in 
our opinion captured by the adiabatic terms linear in dt and V. We will now discuss the 
microscopic basis for Eq. ([9]) containing such terms. 
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V. MICROSCOPIC THEORY OF MAGNETIZATION DYNAMICS 



Once the phenomenological equation for current- driven magnetization dynamics is re- 
duced to the form which requires smooth magnetization variation, slow dynamics, and 
isotropic ferromagnetism, the remaining key questions concern the magnitude and relation 
between the two dimensionless parameters a and f3. The size of the Gilbert damping con- 
stant a is a long-standing open question in solid-state physics, and even a brief review of the 
relevant ideas and literature is beyond the scope of this paper. A recent model calculation 
highlighting the multitude of relevant energy scales that control magnetic damping can be 
found in Ref. Here, we discuss only the ratio (3 /a, since it is of central importance 

for macroscopic current-driven phenomena. As noted above, the ratio (3 /a determines, for 
example, the onset of the ferromagnetic current-driven instability [see Eq. f|T7|) ] as well as 
the Walker threshold current (both diverging when f3/a — 1). The subthreshold current- 
driven domain- wall velocity is proportional to (3 /a [see Eq. fl22l) ]. while (3 /a = 1 is a special 
point, at which the effect of a uniform current density j on the magnetization dynamics is 
eliminated in the frame of reference that moves with velocity v = —V^/sq [see Eq. f llip ]. 
Although the exact ratio is a system-dependent quantity, some quahtative aspects not 
;oo sensitive to the microscopic origin of these parameters have recently been discussed 



13|,|49|,15C 

In Ref. [49|, we developed a self-consistent mean-field approach, in which itinerant elec- 
trons are described by a time-dependent single-particle Hamiltonian 

H = [no + U{r, t)] i + . (H + H.,) (r, t) + n„, (25) 

where the unit matrix 1 and the vector of the Pauli matrices a- = {a^, (^y, o"^) form a basis 
for the Hamiltonian in spin space. Tio is the crystal Hamiltonian including kinetic and 
potential energy. U is the scalar potential consisting of disorder and applied electric-field 
contributions. The total magnetic field consists of the applied, H, and exchange, Hxc, fields. 
Finally, the last term in the Hamiltonian, H^, accounts for spin-dephasing processes, e.g, due 
to quenched magnetic disorder or spin-orbit scattering associated with impurity potentials. 
This last term is responsible for low-frequency dissipative processes affecting a and (3 in the 
collective equation of motion ([9]). 



In time-dependent spin-density-functional theory |44j. 1 771. l78[| of itinerant ferromagnetism. 
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the exchange field Hxc is a functional of the time-dependent spin-density matrix 



p„^(r,t) = (vI/Ur)vI/,(r))f 



(26) 



that should be computed self-consistently from the Schrodinger equation corresponding to 



Focusing on low-energy magnetic fiuctuations that are long range and transverse, we restrict 
our attention to a single parabolic band. Consideration of realistic band structures is pos- 
sible from this starting point. We adopt the adiabatic local-density approximation (ALDA, 
essentially the Stoner model) for the exchange field: 



with direction m = — s/s locked to the time-dependent spin density (1271) (assuming 7 > 0). 

In another simple model of ferromagnetism, the so-called s-d model, conducting s elec- 
trons interact with the exchange field of the d electrons which are assumed to be localized to 
the crystal lattice sites. The (i-orbital electron spins are supposed to account for most of the 
magnetic moment. Because ci-electron shells have large net spins and strong ferromagnetic 
correlations, they are usually treated classically. In a mean-field s-d description, therefore, 
conducting s orbitals are described by the same Hamiltonian fl25|) with an exchange field 
( |28|) . The differences between the Stoner and s-d models for the magnetization dynamics 
are rather minor and subtle. In the ALDA/Stoner model, the exchange potential is (on the 
scale of the magnetization dynamics) instantaneously aligned with the total magnetization. 
In contrast, the direction unit vector m in the s-d model corresponds to the d magnetization, 
which is allowed to be misaligned with the s magnetization, transferring torque between the 
s and d magnetic moments. Since most of the magnetization is carried by the latter, the 
external field H couples mainly to the d spins, while the s spins respond to and follow the 
time-dependent exchange field 0281) . As Axe is usually much larger than the external (includ- 
ing demagnetization and anisotropy) fields that drive collective magnetization dynamics, the 
total magnetic moment will always be very close to m. A more important difference of the 
philosophy behind the two models is the presumed shielding of the d orbitals from exter- 
nal disorder. The reduced coupling with dissipative degrees of freedom would imply that 



Ti. The spin density of conducting electrons is given by 




(27) 




(28) 
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their dynamics are much less damped. (Whether this is actually the case in real systems 
remains to be proven, however.) Consequently, the magnetization damping has to come 
from the disorder experienced by the itinerant s electrons. As in the case of the itinerant 
ferromagnets, the susceptibility has to be calculated self-consistently with the magnetization 



dynamics parametrized by m. For more details on this model, we refer to Refs. [10|, l49 |. 
With the above differences in mind, the following discussion is applicable to both models. In 
order to avoid confusion, we remark that the equilibrium spin density sq introduced earlier 
refers to the total spin density, i.e., d- plus s-electron spin density, while Eq. (1271) refers 
only to the latter. The Stoner model is more appropriate for transition-metal ferromagnets 
because of the strong hybridization between d and s,p electrons. Magnetic semiconductors 
are characterized by deep magnetic impurity states for which the s-d model may be a better 
choice. 

The single-particle itinerant electron response to electric and magnetic fields in Hamil- 
tonian f l25|) is all that is needed to compute the magnetization dynamics microscopically. 
As mentioned above, the distinction between the Stoner and s-d models will appear only 
at the end of the day, when we self-consistently relate m(r, t) to the itinerant electron spin 
response. Before proceeding, we observe that since the constants a and (3 which parametrize 
the magnetic equation of motion (Q affect the linear response to a small transverse applied 
field with respect to a uniform magnetization, we can obtain them by a linear-response 
calculation for the single-domain bulk ferromagnet. The large-scale magnetization texture 
associated with a domain wall does not affect the value of these parameters, in the consid- 
ered limit. The linear response to a small magnetic field is complicated by the presence of 
an electrically-driven applied current, however. Since the Kubo formula based on two-point 
equilibrium Green's functions is insufficient to calculate the response to simultaneous mag- 
netic and electric fields, we chose to pursue a nonequilibrium (Keldysh) Green's function 



formalism in Refs. 



13, 



491]. A technically impressive equilibrium Green's function calcula- 
tion has been carried out in Ref . [s^ , which to a large extent confirmed our results, but also 
contributed some important additions that will be discussed below. 



The central quantity in the kinetic equation approach 



13 



49| is the nonequilibrium 



component of the 2x2 distribution function /k(r,t). In the quasiparticle approximation. 



valid when Axe Ep 



491 ]. the kinetic equation can be reduced to a semiclassical Boltzmann- 



like equation that accounts for electron drift in response to the electric field as well as the 
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spin precession in the magnetic field. The nonequihbrium component of the spin density 
reads s' = {h/2) J rf^fc fk/(27r)3, where fk = Tr[/kCT]: 

9tS —Z X S — Z X U = -- / 777^(Vk ■ <9r)fk . (29) 



h h 2 J (27r)3' K ry K 

s here is the equihbrium spin density of itinerant electrons, Vk = dkB\f_s/h is the momentum- 
dependent group velocity, and the magnetization direction m = z + u is assumed to undergo 
a small precession u relative to the uniform equilibrium direction z. The first term on the 
right-hand side is the spin-current diver gen ce and the last term is the spin-dephasing term 



introduced phenomenologically in Ref. 49|] and studied microscopically in Ref. 13|. As 
detailed in Ref. [43], the spin currents have to be calculated from the full kinetic equation 
and then inserted in Eq. (12^ . The final result (for the Stoner model) is given by Eq. flTU]) 
or, equivalently, Eq. ([9]), with a = p. The latter is proportional to the spin-dephasing rate: 

P = ^. (30) 

The derivation assumes uj,t^^ ^ Axc/^, which is typically the case in real materials suffi- 
ciently below the Curie temperature. The s-d model yields the same result for (3, but 

a = r](3 (31) 

is reduced by the rj = s/so ratio, i.e., the fraction of the itinerant to the total angular 
momentum. [Note that Eq. fl3Tl) is also valid for the Stoner model since then sq = s.] For 
the s-d model, the equation of motion ([9]) clearly cannot be reduced to Eq. (fTOj) . since a ^ [3. 
The steady-state current- driven velocity (122|) for both mean-field models becomes 

f^Vj Vj 

V = = , (32) 

aso s 

where s is the itinerant electron spin density. Interestingly, the velocity fl52]) is completely 
determined by properties of the conducting electrons, even for the s — d model. In the Drude 
model, 

V (x^, (33) 
m* 

where E is the applied electric field, r is the characteristic momentum scattering time, and 
m* is the effective mass of the itinerant bands at the Fermi energy. We expect the velocity 
fl55]) . which is essentially the conducting electron drift velocity, to be suppressed for the 
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s-d model if the d orbitals are coupled to their own dissipative bath, which has not been 
included in the above treatment. 



Ref. 



Sd ] refines these results by relaxing the assumption that Axe ^ Ep and by consid- 



ering also anisotropic spin-dephasing impurities, which results in a ^ j3 for both Stoner and 



s-d models. Ref. 



511] later offered a Keldysh functional-integral approach leading to the same 



results. (These authors also found stochastic torques expressed in terms of thermal fluctu- 



ations ([5]) in the weak current limit; see, however, Ref. [7l| for additional current- induced 
stochastic terms present in the case of an inhomogeneous magnetization.) Consider, for 
example, weak magnetic disorder described by the potential Ha = h(r) ■ & with Gaus- 
sian white-noise correlations {ha{r)hh{r')) oc UaSabS{r — r'), where Ua = U± when a is 
perpendicular (parallel) to the equilibrium magnetization direction. (Spin-orbit interaction 
associated with scalar disorder gives similar results.) For isotropic disorder, U± = U\\, and 
Axe ^ Ep,a/P ~ T] with r] = s/sq, as was already discussed (reducing to 77 = 1 for the 
Stoner model). Even for larger exchange, the correction to this a/ (3 ratio turns out to be 
rather small: For parabolic bands, for example, ~ [1 — {Axc/Ep)'^/A8]r]. This ratio is 
more sensitive to anisotropics U\\ 7^ U±, however, so that in general a/(3 rj even in the 



limit A^jEp ^ |50| 



VI. SUMMARY AND OUTLOOK 

Our microscopic understanding is based on a mean-fleld approximation, in which itiner- 
ant electrons interact self-consistently with a space- and time-dependent exchange fleld. We 
presented results for the local-spin-density approximation and the mean-fleld s-d model. We 
identifled a relation between dissipative terms parametrized by a and (3 and spin-dephasing 
scattering potentials. The central result for the collective low-frequency long-wavelength 
current-driven magnetization dynamics can be formulated as a generalization of the phe- 
nomenological Landau-Lifshitz-Gilbert equation, accounting for the current-driven torques. 
One should in general also include stochastic terms due to thermal fluctuations as well as 



nonequilibrium shot-noise contribution in the presence of applied current j [71j. Despite 
some recent efforts, stochastic effects remain to be relatively unexplored both theoretically 
and experimentally, however. 

The most important parameter that determines the effect of an electric current on the 
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collective magnetization dynamics in extended systems is the ratio (3 /a. We find that 
this ratio is not universal and in general depends on details of the band structure and 
spin-dephasing processes. Nevertheless, simple models give a ~ /5 with the special limit 
a ~ /3 for the Stoner model with weak and isotropic spin-dephasing disorder. Solving the 
magnetization equation of motion for a domain wall is rather straightforward at low dc 
currents, when the wall is only slightly compressed. The domain-wall motion can then be 
modeled within the Walker ansatz, based on parametrizing the magnetic dynamics in terms 
of wall position and spin distortion. Two regimes can then be distinguished: At the lowest 
currents, the wall moves steadily in the presence of a constant uniform current, while above 
the so-called Walker threshold, the magnetization close to the wall center starts oscillating, 
resulting in a singular dependence of the average velocity on the applied current. 

The values of the a and (3 parameters are not affected by the magnetization textures. 
Micromagnetic simulations can provide better understanding of experimental results in the 
regimes where domain walls are not well described by a one-dimensional model. Experiments 
can contribute to the understanding by studying ferromagnets with systematic variations 
of impurity types and concentration, for Py and other different materials. Experimental 
investigation of creep in metalhc ferromagnets at temperatures far below the critical ones, 
as compared to studies [33 



391 ] on magnetic semiconductors close to the Curie transition, 
are highly desirable in order to advance our understanding. 

Besides realistic microscopic evaluations of the key parameters a and (3, the collective 
current-driven magnetization dynamics pose many theoretical challenges, in the spirit of 
classical nonlinear dynamical systems. Current-driven magnetism displays a rich behavior 
well beyond what can be achieved by applied magnetic fields only. At supercritical currents. 



ferromagnetism becomes unstable, possibly leading to chaotic dynamics [79|, although al- 
ternative scenarios have been also suggested SOj. Domain- wall dynamics in a medium with 
disordered pinning potentials pose an interesting yet, at weak applied currents, tractable 
problem. Spin torques and dynamics in sharp walls and the role of strong intrinsic spin-orbit 
coupling (relevant for dilute magnetic semiconductors) are not yet completely understood. 
Oscillatory domain- wall motion under ac currents and in curved geometries is also starting 



35 



8l| . Another direction of 



to attract attention both experimentally and theoretically 

recent activities concern the backaction of a moving domain wall on the charge degrees of 



freedom 



82 



83, 



84 



85 



86|. 
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With the exciting recent and forthcoming experimental developments, the questions con- 
cerning interactions of the collective ferromagnetic order with electric currents will certainly 
challenge theoreticians for many years to come. The prospects of using purely electric means 
to efficiently manipulate magnetic dynamics are also promising for practical applications. 
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